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Two improved superposition schemes for the construction of approximate multi-caloron- 
anticaloron configurations, using exact single (anti)caloron gauge fields as underlying building blocks, 
are introduced in this paper. The first improvement deals with possible monopole-Dirac string inter- 
. actions between different calorons with non-trivial holonomy. The second one, based on the ADHM 

formalism, improves the ( ant i-) self duality in the case of small caloron separations. It conforms with 
Shuryak's well-known ratio-ansatz when applied to instantons. Both superposition techniques pro- 
vide a higher degree of (anti-)selfduality than the widely used sum-ansatz, which simply adds the 
(anti)caloron vector potentials in an appropriate gauge. Furthermore, the improved configurations 
(when discretized onto a lattice) are characterized by a higher stability when they are exposed to 
lattice cooling techniques. 
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I. INTRODUCTION 

o ■ 

Lately, calorons with non-trivial asymptotic holonomy [l|, 0, H[ (besides the longer known ones 
with trivial asymptotic holonomy, HS calorons have attracted some interest, both as solutions 
themselves Q and as underlying building blocks of some semiclassical-like modeling of QCD in the 
finite temperature case Q by a caloron-anticaloron gas model. The motivation of such simulations 
was to check whether these solutions, rather than instantons, are able to produce confinement. In 
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some sense these solutions are better adapted to varying physical situations including the confinement 
phase where the average holonomy is non-trivial. One interesting feature of these solutions is that 
they can occur decomposed into constituents with fractional topological charge. As we will see, this 
fact causes part of the problems that require to reconsider the superposition of calorons. 

The semiclassical picture of QCD dates back to Callan, Dashen and Gross (CDG) [8|. The authors 
have worked out the instanton effects on the heavy quark potential @ and proposed a confinement 
mechanism for quarks being enclosed in bubbles ( "bags" ) of the wrong, perturbative vacuum [l(| [H| • 
Finally, this mechanism could not be verified by later research. It has become also an accepted 
fact, that the instanton gas or liquid is unable to generate a confining potential between static color 
charges. The original proposal of CDG, however, to explain confinement by merons [l2j is fascinating 
people working on confinement until today [13]]. Merons are sort of singular half-instantons eventually 
liberated by a phase transition from weak to strong coupling. Each of them is carrying half-integer 
topological charge. Since then, the concept of fractional topological charges ("instanton quarks" 
of charge Q = l/N co i OI ) is vivid in the literature (see Refs. [3, and Ref. [l6| and references 
therein). Some indications for the existence of smooth Q = 1/2 lumps in SU(Z) lattice configurations 
and/or their relation to confinement have been published in the 90-s [l?], EH, Gil H3, HH, [HJ and more 
recently (23| . The corresponding solutions for zero temperature are still elusive. 

The application of instanton physics to finite temperature QCD was mainly influenced by the 
seminal paper [241 ]. Therein, for some reasons, the potential significance of Euclidean solutions with 
non-trivial holonomy has been overlooked. Therefore the instanton model (for a general review see 
Ref. J25|) was applied to T ^ exclusively based on HS calorons [26|. Diakonov has recently pointed 
out [27J that these arguments are not compelling for dense multi-caloron-anticaloron systems. His 
paper has stimulated more attention to these new caloron solutions. Ideas about caloron constituents 
as relevant degrees of freedom for the confined phase from the point of view of electric-magnetic duality 
have been developed in Ref. [28[ . 

The encouragement for the simulation of gases (or liquids) of calorons and anticalorons, apart from 
the motivation described above, came from lattice studies. Soon after the discovery of the caloron 
solutions Garcia Perez et al. p9l [30j have considered such objects on a finite lattice, partly with 
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twisted boundary conditions. Later on, in Refs. [3l|, [12, [H, [Hj] the possibilities to detect calorons in 
generic Monte Carlo configurations have been developed further, in a similar spirit as earlier instanton 
(caloron) searches at zero (non-zero) temperature [35l [36l . [37| . In view of this development it was 
natural to develop the tools needed to model multicaloron systems. The necessity of simulating such 
systems reflects the obvious fact that a strictly semiclassical calculation, which has been done for 
a single caloron (38|, is not feasible in a strongly coupled system like QCD near to the deconfining 
temperature. 

The mathematical-physical interest in the solutions themselves has led to remarkable success p| @ in 
recent years concerning the construction and description of selfdual gauge fields with higher topological 
charge Q. The parametrizations for \Q\ ^ 1 are nevertheless restricted to cases of special symmetry. 
A general parametrization is still lacking. Apart from this, the modeling of gauge field trajectories 
of QCD in a semiclassical-like fashion in terms of calorons is not using (like earlier instanton liquid 
models [39|, 0, |4l[) exact solutions of arbitrary topological charge but approximate, only piecewise 
perfect solutions. The approximate solutions should be superpositions of elementary building blocks 
with Q = ±1. Of course, there is necessity and room to perfectionate the superposition scheme such 
that the action and the topological charge of purely selfdual superpositions would be additive to a 
high precision. This is not the case for the simple sum-ansatz in terms of the vector potentials A$ (x) 
of single calorons and/or anticalorons, except for an infinitely dilute gas. For the case of calorons, this 
improvement is addressed here for the first time. 

The first problem of any superposition of calorons with definite holonomy leads us beyond the simple 
sum-ansatz in terms of the vector potentials A^\x) (with i labelling the individual (anti)calorons). 
There is the so-called algebraic gauge where the additive superposition is indeed possible. The gauge 
transformation to the periodic gauge, which is the more appropriate one to build the ensemble of 
model configurations, requires all of them to have the same asymptotic holonomy [5j]. Once this is 
adopted, a reasonable model of QCD can be constructed 0]. 

We introduce in the present paper two new techniques to facilitate the further improvement of such 
superpositions. One of them is specific for calorons with non-trivial holonomy, which may exist in 
a state "dissociated" into monopole constituents. These monopoles are connected by non-Abelian 
strings of strongly peaked, but regular vector potential with low action density. They turn into 
usual Dirac strings, being Abelian and singular, when the monopole separation tends to infinity. For 
simplicity we will denote the regular, non-Abelian string as "Dirac string". The presence of such 
"Dirac strings" creates problems whenever one attempts to add another caloron [f| . Since the number 
of constituent monopoles in a caloron is equal to the number of colors, N co \ or , the case of SU(2) 
gluodynamics with two constituents and one Dirac string per caloron is a convenient testing ground 
for this type of problems. 

The other improvement problem is connected with caloron interactions arising when constituents 
(lumps of action) happen to come too close to each other when the collective coordinates are sampled. 
For these cases, the ADHM formalism (42| offers the necessary intuition for the improvement of the 
superposition. 

We should remark that the quest of improvement is restricted only to semiclassical-like simulations. 
Sampling configurations with merons and instantons in the regular gauge [l3| or with more general 
pseudoparticle fields [43] do not need to approximate these by solutions of the Euclidean equations 
of motion. Superpositions in this case have to be subjected to importance sampling which induces 
strong correlations. It is these long range correlations which generates confinement in this class of 
models. Our approach is complementary in the sense that our superpositions are much closer to real 
solutions and can be randomly sampled in moduli space Q- The long range correlation required for 
confinement in this case is due to the non-trivial holonomy tail of the gauge potentials. 

This paper is organized as follows. In Section [TT] we briefly recall the properties of single 
(anti)calorons. In Section IIIII we explain under which circumstances, in the algebraic gauge, an 
additive superposition is possible. There we also introduce technical details like the use of the lattice 
discretization and some observables characterizing the quality of the constructed superposition. We 
explain e-improved cooling for later use in this paper. In Section HVl we propose a solution for the ob- 
struction problem presented by Dirac strings and demonstrate the effect for a particular two-caloron 
superposition. In Section [V] we describe the ADHM-inspired improvement for closely overlapping 
caloron constituents. The success is demonstrated for a particular four-caloron superposition. In 
Section IVII we summarize and draw conclusions for further work. 
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II. CLASSICAL SOLUTIONS OF SU{2) YANG-MILLS THEORY 

The instanton, discovered in 1975 [44]], is a classical solution of the Euclidean Yang-Mills equation 
of motion in flat 1R 4 , representing zero temperature, with localized action density and one unit of 
topological charge. For any number of colors, the instanton is basically an SU(2) object that is 
parametrized by 4 • N co \ OT collective coordinates. For iV co i or = 2 these are its four-dimensional 
position in space-time, a size parameter and three parameters describing a global rotation in SU(2) 
group space. For higher N co \ m the additional parameters describe the embedding of SU(2) into 
SU(N co \ OI ). Classical solutions with higher topological charge can in principle be constructed by 
means of the ADHM formalism [42j], although they are not analytically available in general. 

The analog of the instanton in the case of non-zero temperature T = 1/(3 ^ was discovered soon 
after the T = instanton by Harrington and Shepard in 1978 0] and was accepted as an appropriate 
semiclassical background at finite temperature. It will be referred to as 'HS caloron' throughout the 
present paper. 

Twenty years later, Kraan and van Baal P, and Lee and Lu [|| extended the parameter space of 
the HS caloron by an additional parameter, the asymptotic holonomy Vooi which is defined as 

7>co - e 2 ^ ? = Jim P(x) , (1) 

\x\ — >oo 

where f denotes the vector of Pauli matrices and 

P{x) = P exp J A 4 {x,t)dt^j G 577(2) , (2) 

where path ordering is implied by the P symbol. According to this definition the Polyakov loop 
£(x) = ^TrP(x) at spatial infinity is directly determined by the parameter uj. The generalized 
solution, which we will call KvBLL caloron, is reduced to the HS caloron for trivial asymptotic 
holonomy, taking values in the center of the gauge group, Voo G Z{2). In the case of general iV co i or , 
this solution can also be constructed ^|. Generically, it is not a simple embedding of an SU{2) 
solution into the bigger group. The present paper, however, is restricted to the case of N co \ or = 2. 

The KvBLL caloron with arbitrary asymptotic holonomy is a selfdual gauge field, i.e. = F^ v 
with Fn„ — ^€t_tva0F a f}, solving the Euclidean equation of motion at finite temperature. It was first 
constructed in the so-called algebraic gauge by means of the ADHM formalism [42| by Kraan and van 
Baal [3, [2] and by Lee and Lu [3] as a selfdual gauge field in flat K 4 with the periodicity condition 
A°j 9 (x+/3it) = VooA^ 9 (x)T > l , where it denotes the unit vector along the time direction. The resulting 
vector potential A°j 9 defined on 1R 4 can then be transformed into a selfdual and time-periodic gauge 
field Al er , defined on R 3 x S 1 , by a gauge transformation f2(x) non-periodic in time, according to 

A£ er (x,imod/3) = Q,{x)A^ 9 {x,t)n\x)-iQ,{x)d^{x) , (3) 

with x = (x, i) and 

n(ar) = e' 2 "™* 7 ' 9 . (4) 

This is the KvBLL caloron in the periodic gauge with arbitrary asymptotic holonomy 
"Poo = exp(27ricJr) G SU(2) and topological charge Q = 1. One finds that the KvBLL caloron is 
also described by 8 collective coordinates in addition to its holonomy parameter uj. However, the 
physical interpretation of these collective coordinates is different from that of the instanton parame- 
ters, as will be discussed in the following. 

One can imagine the KvBLL caloron in the algebraic gauge as an infinite chain of equally separated, 
identical instantons in flat 1R , aligned along the time direction with each subsequent instanton rotated 
by the holonomy Voo relative to the preceding one - as long as the size p of the instantons is small 
compared to the inverse temperature (3. In that case it is appealing to parametrize the caloron in 
terms of the collective coordinates of the instanton. For p <C /3 the corresponding caloron consists of 



FIG. 1: Action density (left) and Polyakov loop field £(x) (right) of a KvBLL caloron with maximally non- 
trivial holonomy \uj\ — 0.25, p = 1.0 and [5=1, where the action density is given in instanton units Si. For 
sufficiently large p/[3 the two constituents of the KvBLL caloron become well separated. Here, pzD denotes 
the radius of the almost 0(3)-symmetric monopole action lumps in 3 dimensions. 

one lump of action, which is approximately 0(4)-rotationally symmetric in space-time, and its radius 
can be described by the parameter p. 

In the opposite case, p ~ j3 or p 3> /3, the caloron is dissociated into two constituents, as shown in 
the action density plot in Figure [TJ For p — > oo the action density of these constituents becomes static 
in time, breaking down the prior 0(4) symmetry to an 0(3)-rotational symmetry in 3-dimcnsional 

(1 2) 

position space, and the corresponding 3-dimcnsional radii p 3 p of the two constituents converge to con- 
stants approximately given by p^ ~ /3/(4oj) and p^ w /?/ (4cD), respectively, where to = \uj\ € [0; 1/2] 
and Q= 1/2 — u>. The constituents can be identified with BPS monopoles [46[ each carrying the 
fraction 2lo or 2cD, respectively, of the total action and topological charge. The distance between the 
monopole positions zi, z 2 is given by d = \zi — z 2 \ — irp 2 / ' [3. Due to the emergence of these two 
monopoles at larger p it seems more natural to use their positions Z\ and z 2 for the parametriza- 
tion of the caloron instead of the instanton parameters involved in the construction scheme. The 
corresponding vector potentials in the algebraic and the periodic gauge are analytically given in [H, • 
Three further remarks shall be given here. 

(i) In the case of non-trivial asymptotic holonomy Voo the caloron gauge field does not vanish at 
spatial infinity. In the periodic gauge the A± component of the vector potential converges to a constant 
at spatial infinity whereas the remaining three components approach zero, 

Jim Af» = 2™r//3, lim A^ 3 (x) = . (5) 

\x\ — >oo \x\ — >oo 

The non- vanishing vector potential will become an issue when one tries to superpose KvBLL calorons. 

(ii) In the case u> — > or uj — > the asymptotic holonomy becomes trivial, one constituent vanishes 
and the remaining one becomes a HS caloron. 

(Hi) Finally, the anticaloron, which is the anti-selfdual analogue of the caloron, carrying negative 
topological charge, can be obtained from the caloron gauge field by 

A? nt *(x,t) = -Ai(-x,t), 1=1,2,3, 

AT u {x,t) = A 4 (-x,t) . (6) 

III. SUPERPOSING CALORONS: A FIRST APPROACH AND ITS PROBLEMS 

In general, the superposition of N sclfdual gauge fields A^\x), i — 1,...,N, does not provide a 
selfdual gauge field configuration due to the non-linearity of the field strength (x) in terms of the 
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vector potential. However, since the non-linearity arises from the commutator term [A^{x), A v (x)\ 
appearing in the expression for the field strength, such a superposition has the potential of approxi- 
mately describing a selfdual gauge field, provided that at every space-time position x the superposition 
is dominated by the contribution A^} (x) of one single caloron, or that all vector potentials are almost 
Abelian. 

This well-known observation has been exploited in numerous semiclassical model calculations based 
on instantons in the singular gauge Q relying on the strong localization of the instanton gauge field 
in space-time. For instantons the simplest possible superposition scheme, which plainly consists of 
adding all involved vector potentials according to the simple sum-ansatz 



can successfully be applied as long as the instanton cores do not come too close to each other. The 
meaning of the sloppy formulation 'not too close' was clarified in [47] . establishing the quantitative 
criterion p\p\j\t,\ — ^2 1 4 < 1/200 for the applicability of the simple sum-ansatz for the superposition 
of two instantons with sizes p\, P2 and positions £1, £2- 

When dealing with calorons with non-trivial holonomy, however, this criterion no longer holds, 
since the KvBLL caloron gauge field then does not vanish even far away from the monopole locations 
according to Eq. ([5]). 

Applying naively the simple sum-ansatz ([7]) would induce an interaction between calorons even at 
infinite separations. This is the first caloron-specific problem appearing in superpositions. It can 
easily be overcome by superposing the caloron vector potentials in the algebraic gauge, in which the 
potential A^ 9 vanishes at spatial infinity such that the sum-ansatz does no longer yield interactions 
between infinitely far separated calorons Q. On the other hand, the resulting gauge field is then no 
longer periodic. Periodicity has to be restored by applying a specific gauge transformation to the sum 
of the vector potentials. This leads to the modified ansatz 



which we will refer to as sum-ansatzin the following, where (3 = 1 was chosen for convenience. However, 
this superposition scheme is only valid under the condition that only calorons and anticalorons with 
identical holonomy be superposed because the gauge transformation necessary to bring the caloron 
vector potential from the algebraic to the periodic gauge depends on the holonomy according to 
Eqs. (0 2]). It would not be possible to restore the periodicity of a configuration of calorons first 
constructed in the algebraic gauge with a global gauge transformation, if the configuration were a 
superposition of calorons with different holonomies [f|. 

An additional caloron-specific issue concerning the superposition of well-dissociated calorons with 
non-trivial holonomy, which cannot be dealt with by means of the sum-ansatz ([5]), is due to the 
existence of the "Dirac string" , which in our context is a string of strongly peaking, but regular vector 
potential connecting the two monopoles within a caloron. Despite its strength in A^, the Dirac string 
is almost free of action thanks to a fine-tuning of its gauge field. However, there can potentially 
be strong interactions between this Dirac string and other calorons or other monopoles, respectively, 
spoiling the selfduality of the superposed configuration. A more appropriate superposition scheme 
taking care of the caloron-Dirac string interactions will be presented in Section IIVI Furthermore, 
similar to the case of instantons, the sum-ansatz also fails, if two calorons or monopoles come too 
close to each other. An improved superposition scheme for the latter case, leading to smaller deviations 
from selfduality, is introduced in Section [VI 

So far, our discussion was referring completely to the continuum. However, for the discussion of 
the improved superposition schemes to be presented we shall find it quite instructive to apply certain 
lattice techniques to the gauge field configurations we have generated. For this purpose, we introduce 
a grid (called "lattice" in the following) and represent the gauge field A^(x) of the caloron ensemble 
by a set of link variables according to 




(7) 




(8) 




(9) 



G 



where the path ordering symbol P will approximately be accounted for by dividing each path of 
integration into n sub- intervals with their lengths a/n -C p being small compared to the average 
caloron size, such that the vector potential inside these sub-intervals can be considered to be constant. 
Furthermore, the convergence of the link variables is ensured by adjusting the number of sub-intervals 
dynamically, corresponding to the variability of the gauge field, and checked by comparing the resulting 
parallel transporter for different values of n. 

This introduction of a lattice makes the constructed gauge field configurations accessible to cooling 
procedures [H, E§| ■ From the cooling history we can then assess the stability of the configuration 
under the cooling process, which will give information on the "distance" to the "nearest" classical 
solution in the space of all field configurations, and thus helps to evaluate the quality of the applied 
superposition scheme. Here, we will use the so-called e-cooling [23|, [5(| with e = 0, based on an 
improved lattice action <Si mpr [A] = S e =o[-<4], the lattice artefacts of which only start in 0(a 4 ). 

To estimate the importance of a constructed configuration according to its weight contributing to 
the partition function Z it is interesting to consider its total action. We define the global quantity 
7 as the ratio of the total lattice action S[A] (in instanton units Si = Sir 2 / g 2 ) and the topological 
charge Q[A] according to 



7 \Q[A}\> 

where we also use an improved definition of the lattice topological charge Q[A] [5l| . 

Finally, as a measure of the local violation of (anti-)selfduality we introduce the quantity 



(10) 



S(x) 



E 



FyL V (x) - K F fiv (x) 



> 



with 



sign(Q) = ±1. 



(11) 



In what follows 8{x) will directly be determined from the continuum gauge fields. All other results 
presented for the action, topological charge, and their densities will be measured by the lattice tech- 
niques introduced above. Both the quantities 7 and 8 (x) reflect the deviation from (anti-)selfduality 
for k = +1 (— f). For purely (anti-)selfdual configurations A^x) the local measure S(x) is identical 
to zero, while the global measure 7 becomes unity. Of course, these definitions make only sense for 
superpositions built exclusively out of either calorons or anticalorons in an attempt to construct higher 
charge calorons or anticalorons. 



IV. REDUCING THE CALORON-DIRAC STRING INTERACTION 



In this Section we deal with the interaction between the Dirac string of a dissociated caloron 
and any further (anti)calorons. Due to its strength the Dirac string has the potential of creating a 
tremendous amount of unwanted additional action and of completely spoiling the (local) selfduality of 
a multi-caloron system created by the sum-ansatz. Certainly, this problem becomes most prominent 
if a further monopole or caloron, respectively, is placed between the two monopoles of a dissociated 
caloron, i.e. on top of the Dirac string (see also |5|). 

We will see that this problem can be cured, at least if the constituents causing the Dirac string are 
sufficiently far separated from each other, by applying a certain rotation in color space (not a gauge 
transformation) on the caloron located on top of the Dirac string. This result will be obtained from the 
observation that the Dirac string of a caloron with non-trivial holonomy is almost pure gauge for well- 
dissociated calorons, i.e. the vector potential of the Dirac string between the monopoles can almost 
be gauged to zero. However, the erased components of the vector potential will reappear outside of 
the caloron. Finally, the improved superposition scheme will be compared with the sum-ansatz. 



A. Structure of the Dirac string vector potential 



We begin with studying the structure of the Dirac string gauge field for a single, well-dissociated 
caloron, i.e. p 3> /3, in the algebraic gauge and with maximally non-trivial holonomy to = 0.25. For 
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convenience we fix u)f = ujt 3 and align the orientation of the monopole axis along the a^-axis with 
the caloron center position located at the origin. We now consider the three Lie algebra components 
of the caloron vector potential 



A^x) = A%\x)- 



(12) 



at positions x — (x, t) close to the connecting axis and between but not too close to the two monopole 
positions. In the large p limit one easily derives from the analytical equations of the caloron gauge 
fields (given in [l], Q), that the vector potential A^{x) is dominantly Abelian and governed by the 

(31 

A\ 2(5) components according to 



A3 Ax) 



< 




< 





(13) 
(14) 



In fact, one can show that the A^l(x) components are proportional to 1/(A + e) in the large p 

limit, where A = \/ x\ + x\ is the distance to the monopole axis and e = e(x,p) is a small number, 
exponentially suppressed by the p parameter. This is what we referred to with the term "strongly 
peaking, but regular" in the preceding sections. Furthermore, the gauge fields become static in time 
and in the region between the monopolcs also independent of the a^-coordinate 



A[ 3) 2 (x) « A\%( Xl ,x 2 ,x 3 ,t = 0) 



,(3) 



Af\{i) « A^{ Xl ,x 2 ,x 3 = 0,t) 



.(3) 



(15) 
(16) 



These observations are illustrated in Figure [5] showing for the situation under discussion the behavior 
of the three Lie algebra components of the corresponding caloron vector potential in the xi-X3-plane. 

We will now find a gauge transformation that puts the dominant components A\ 2 (x) between the 
constituents to almost zero. From the above observations we conclude that the transformation 



A Ax) = e iG ^ ■ A^x) ■ e - iG ^ - d^G(x)r 3 , 



(17) 



which we refer to as DS transformation, with a function G(x) = G(x\,X2) depending only on the x\- 
and a^-coordinate, will do the job provided that G{x\, x%) fulfills the two conditions 



d 1 G(x 1 ,x 2 ) = A\' (x 1 ,x 2 ,x 3 = 0,t = 0) 
d 2 G(xi, x 2 ) = A\ '{xi,X2 i x 3 = Q,t = 0) 

1 (3) 



(18) 
(19) 



This guarantees that at x 3 = the formerly dominating A\ 2 components are completely erased. 
However, the existence of such a function G is non-trivial. The above requirements Eqs. (|18j) and 
(fl9| have to be consistent with the condition of integrability d\d 2 G(x\, x 2 ) = d 2 diG(xi,x 2 ). A first 
indication of the existence of the desired function G can be obtained from a simple consideration 

taking the vanishing action density of the Dirac string and its Abelian dominance into account. 

(3) 

However, an explicit calculation of the vector potentials A\ ^(x) at x 3 = and t = in the large p 
limit for A = \J x\ + x\ <C T 



■ z 2 \, where z\ and z 2 denote the monopole positions, 



clarifies this question. It yields 



A ( f ) (x 1 ,x 2 ,x 3 =0,t = Q) 
A 2 3) (x u x 2} x 3 = 0,t = 0) 



_X2 (I _ U 
A \X 2r 2 
xi (1 _ 1 A_ 
A I A 2 r 2 



O 



o 



(20) 
(21) 



While the first summand is integrable as presumed, the second summand is not. However, it is sup- 
pressed compared to the first one for small distances A to the monopole axis. Using polar coordinates 
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FIG. 2: The three Lie algebra components of the Euclidean vector potential A^(x) in the algebraic gauge 
are shown for a dissociated caloron with maximally non-trivial holonomy |w| = 0.25, p = 2 and (3 = 1. The 
monopoles are located on the a;3-axis. 



{x\,X2) — > (A, a) to parametrize the X\-X2 plane at 2:3 = the requested function G can now be 
denned sufficiently close to the monopole axis by the path integral 

G(X,a) = r-^ 8 ^? V f^^'"''"^ '^^ V=°°« + ^r^ • (22) 

J V Acos(a') ; yA^(X,a',x 3 =0,t = 0) J \r 2 J K ' 

where adequate mappings between polar coordinates (A, a) and Cartesian coordinates (£1,2:2) are 
implicit. One sees that this function G is not continuous. Close to the axis the function value has 
a 2tt discontinuity at a = 0. The expression e lGT3 , however, is continuous for sufficiently small A. 
Hence, the Dirac string vector potential for < A <C r, s is almost pure gauge in the limit p — > 00. 

So far, we were only referring to the case of maximally non-trivial holonomy, i.e. \uj\ = 0.25. For 
general non-trivial holonomy \uj\ ^ 0.25 the ratio of the 3-dimcnsional monopole radii is given by 
cj/u), such that for a fixed value of p the constituents start to overlap, as soon as lo or u), respectively, 
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becomes too small. In that case the DS-superposition scheme, to be presented in the following, is no 
longer reasonable. However, even for small values of uj or Q the constituents become well-separated 
again, if the value of p is sufficiently increased. In an adequately chosen region between the two 
constituents the vector potential then again obeys the relations (|131 [T4l [T51 ITB|) . This can easily be 
seen from the so-called far-field limit 

Affix) = -\^l v d v ln^/(z) with tff{x) = + ^ (23) 

which describes the caloron gauge fields at sufficiently large distances away from the monopole cen- 
ters [5(. Here, the same assumptions and notations as in Section [TTJ are used and uj = The 
far-field is completely Abelian and in the algebraic gauge it is independent of uj and hence, one can 
easily derive from (|23|) the same typical structure of the Abelian part of the Dirac string that was 
given in (|20|21|1 . The DS-superposition scheme, to be introduced in the following, is therefore also ap- 
plicable (with the same function G{x\,x>2)) for general non-trivial holonomy with the only restriction 
that for a given value of the parameter p the region between the constituents where the vector poten- 
tial behaves as described in (|20|2ip is decreased the more the holonomy deviates from its maximally 
non-trivial setting. 



B. Application to caloron superpositions 



We now introduce an arbitrary, but localized object A 0h: > [x) to be located on top of the Dirac string 
of the previously considered caloron A Cal (x), which is still assumed to be centered at the origin with 
the monopole axis aligned along the ^-direction for convenience. The superposition scheme resulting 
in the improved gauge field Az (x) for the superposition of both gauge fields is now worked out in 
three steps. 

(i) The DS gauge transformation Eq. (fT7|) is applied to A Cal (x) according to 



A'f al (x) = e tG{x)T * ■ A Ca \x) ■ e ~ lG{x)T3 - d^G(x) 



?3 



(24) 



where G(x) = G{x\, £2). This transformation reduces the vector potential of the Dirac string in the 
space between the two constituents. At the same time it creates two new Dirac strings of infinite 
length outside of the caloron which begin at the monopole locations as shown in Figure [3l 



E I4.C»0I 

10 r 




FIG. 3: The square norm ^ I^Ai( a; )| 2 01 the vector potential in the algebraic gauge before (left side) and after 
(right side) applying the DS gauge transformation Eq. (|17[) to the strongly dissociated caloron A^ al (x). The 
monopoles of the latter are located on the a?3-axis. By the transformation the gauge field in the space between 
the monopoles is approximately put to zero, while two new Dirac strings are created outside the caloron again 
aligned parallel to the x^-axis. 
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(ii) The gauge fields, both brought to the algebraic gauge, are added resulting in the interim result 
j4„(x) given as 

A , li (x)=A'™(x)+AW(x). (25) 

(Hi) The inverse gauge transformation is applied to A^x) to remove the newly created Dirac strings 
and put them back in place, i.e. back into the space between the two monopoles, yielding the improved 
gauge field A® (x) of this DS- superposition scheme according to 

A° s (x) = e - iG(x)T3 ■ A'^x) ■ e lG{x)T3 + d a G(x)r 3 

= A Cal (x) +e - lG ^ ■ A° b ^(x) ■ e iG ( x)Ts , (26) 

unchanged rotated 

where the appropriate gauge transformation to the periodic gauge is implicit. 

To summarize, the three steps (i) - (Hi) can all be combined together by simply gauge rotating 
(not gauge transforming) the vector potential A Gb: > (x) of the object to be located on top of the Dirac 
string before adding it to the unchanged field A Cal (x) . 

Three further remarks are in order here. Firstly, this ansatz is also valid for anticalorons. Due 
to the transformation prescription Eq. ((6|) for the construction of anticalorons out of calorons the 
appropriate function G dealing with the Dirac string of the anticaloron even is the same. Secondly, 
this technique can be generalized to an arbitrary number of Dirac strings with arbitrary positions 
and spatial orientations. For this purpose it is sufficient to apply an adequate number of such gauge 
rotations (one for each Dirac string to be dealt with) on the involved monopoles/calorons. And finally, 
the latter statement is true even if situations arise in which the two monopoles of a well separated 
caloron would need to be gauge rotated differently. Such a situation arises, for example, in the case 
of two intertwined calorons as sketched in Figured] for A/i = 0. In this case the two inner monopoles 
would need to be rotated while the two outer ones must not be gauge rotated. However, we have seen 
that the Dirac string connecting the two monopoles consists only of the lot component. Since the 
needed gauge rotation is exactly around this axis in colour space, the rotation has no effect on the 
Dirac string. Hence, the rotation can effectively be stopped at any point of the Dirac string, allowing 
the two monopoles to be rotated independently of each other. Exactly this will be demonstrated in 
the next section. 



C. Comparing DS-improved solutions with the sum-ansatz 

The DS superposition scheme Eq. (|26p introduced above is now compared with the sum-ansatz (|8|). 
Both gauge fields A® s (x) and A s ^ m (x), so far given in analytical form, will be discretized on a lattice 
by a set of link variables according to Eq. (|9|), thus making the gauge fields accessible to cooling 
procedures. The (in) stability of the improved and the unimproved configurations during the cooling 
process will give information about how close the obtained gauge fields are to classical configurations. 

For the comparison a simple composition made out of two intertwined calorons (Ah = 0) as sketched 
in Figure[2]is used. Both calorons are well-dissociated (p — 1.5, (3 — 1), and each of them penetrates 
with its Dirac string one monopole of the other caloron. Certainly, one expects the selfduality to be 
completely destroyed at the inner monopole positions. 

Figure [5] shows the action density, the topological charge density, the Polyakov loop field £(x) and 
the local violation of selfduality S(x), as defined in Eq. (fTTj) , for the configuration obtained from the 
sum-ansatz. The equation of motion is strongly violated at the locations of the inner monopoles as 
expected. This can be seen looking at the local violation of selfduality 5(x) as well as deduced from the 
fact that topological charge density and action density do not match at all in these regions. Moreover, 
the violation is strong enough to have strong impact also on the global measure 7, i.e. the ratio of 
overall action and topological charge as defined in Eq. (fT0|) . which should be unity for exactly self dual 
configurations. From Table U we find 7 = 5.7 ^> 1. 
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FIG. 4: Sketch of the caloron composition to be used for the comparison of the superposition schemes. The 
two calorons are well-dissociated (p = 1.5, (3 = 1), and carry maximally non-trivial holonomy | cZ3 1 = 0.25. The 
distance between their (parallel) monopole axes is denoted as Ah. We will consider the case Ah = 0, where 
one obtains a pair of intertwined calorons penetrating each other with the Dirac string, and Ah ^ 0. The signs 
(±) correspond to the magnetic charge of the monopoles and the straight lines symbolize the Dirac strings. 



Action density Topological charge density 




FIG. 5: Action density in instanton units, topological charge density, Polyakov loop field l{x) and the local 
violation of selfduality 5(x) for the gauge field configuration obtained by the sum-ansatz from the caloron 
composition with Ah = as sketched in Figure [4] 

In order to reduce the violation of selfduality e-cooling with e = is applied to this gauge field 
configuration. Table |T] shows the corresponding cooling history, i. e. the number of cooling steps 
together with the corresponding values for the action, the topological charge and 7. By the strong 
change of the values for the total action and topological charge one sees that the cooling process leads 
to essential alterations of the configuration originally provided by the sum-ansatz. 

In Figure [S] the same composition of calorons is shown but this time after the DS improvement has 
been applied. One sees that now the action density and the topological charge density do approxi- 
mately match and also the local violation of selfduality S(x) is by a factor O(10 3 ) smaller compared 
to the unimproved case. Furthermore, we learn from Table [TT] that the overall action surplus as ob- 
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Cooling steps 


Action S[A] 


Topological charge Q[A] 


1 = S/Q 





9.209 


1.612 


5.714 


1 


5.046 


1.643 


3.072 


5 


2.710 


1.836 


1.476 


10 


2.124 


1.934 


1.098 


20 


1.972 


1.935 


1.020 


50 


1.938 


1.933 


1.002 


100 


1.934 


1.933 


1.001 



TABLE I: Cooling history of the gauge field configuration obtained by the sum-ansatz from the two-caloron 
composition with Ah — as sketched in Figure |U The total action S[A], the total topological charge Q[A] as 
well as the ratio 7 = S[j4]/Q[A] are shown after certain numbers of e-cooling steps with e = 0. 

tained from the quantity 7 is approximately 6 %, thus much smaller than for the previously discussed 
sum-ansatz. Besides this observation we can conclude from the cooling history presented in Table HT1 
that the DS-improved configuration already is relatively close to a classical configuration, since there 
occurs only little alteration of the configuration during the fast convergence of the cooling process. 




FIG. 6: The same as Figure [5] for the gauge field configuration obtained by the DS-improved superposition 
scheme from the caloron composition with Ah = as sketched in Figure [4] 

Finally, let us answer the question how much the benefit of the DS-improvement depends on the 
transverse separation Ah between the monopole axes of the two calorons. For large separation Ah — > 
00 both superposition schemes should be equivalent. In Figure [7] we show the dependence of the total 
action and topological charge on the separation Ah for both superposition schemes. As expected, the 
sum-ansatz and the DS-improvement yield the same results for large separations Ah. For Ah < 2/3, 
which corresponds to the case when Dirac string and monopole touch each other, since the 2D-radius 
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Cooling steps 


Action S[A] 


Topological charge Qb4 M ] 


1 = S/Q 





2.073 


1.941 


1.068 


1 


1.969 


1.940 


1.015 


5 


1.940 


1.937 


1.001 


10 


1.936 


1.936 


1.000 


20 


1.934 


1.934 


1.000 


50 


1.933 


1.933 


1.000 


100 


1.931 


1.932 


0.999 



TABLE II: Cooling history with e = for the gauge field configuration obtained by the DS-improved su- 
perposition scheme from the caloron composition with Ah = as sketched in Figure [4] to be compared at 
corresponding step number with Table [I] 



of the first and the 3D-radius of the latter are approximately /3, the sum-ansatz produces strong 
Dirac string-monopole interactions and thus a great action surplus, while the DS-improved results are 
almost unaffected by the transverse displacement. The latter observation also applies to the (lattice) 
topological charge of the improved configuration. The charge differs by « 5 % from Q = 2 even for 
asymptotic displacement which is a finite- volume or discretization effect. For the unimproved case the 
topological charge deviates much stronger from Q = 2 for Ah < f3, which is due to the finite lattice 
discretization and the vector potential becoming rather unsmooth in the interacting region. 

We can conclude that the major part of the action surplus induced by caloron-Dirac string inter- 
actions can successfully be removed by the DS-superposition scheme, which becomes important at 
monopole-Dirac string separations of Ah — 2(3. Alternatively, we have also seen that this action 
surplus can also be reduced by applying a moderate number of cooling steps. Due to the tremen- 
dous strength of the caloron-Dirac string interaction such a cooling procedure would definitely be 
a necessary step in the semiclassical modeling based on well-dissociated calorons superposed by the 
sum-ansatz. However, whereas the application of a moderate number of cooling steps to a multi- 
caloron system seems to be acceptable, as long as the caloron or monopole radii are large compared 
to the lattice spacing, it becomes quite problematic for smaller calorons, since these small objects 
tend to fall through the grid during the cooling process. Cooling is similarly problematic for mixed 
caloron-anticaloron superpositions. The DS-superposition scheme therefore provides a very interest- 
ing technique for the construction of multi-caloron systems containing strongly dissociated calorons 
as well as small calorons, both of selfdual and anti-selfdual character. 
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FIG. 7: Dependence of the total action S[A] in instanton units Si (left) and of the total topological charge 
Q[A] on the transverse distance Ah between the monopole axes for the two-caloron composition sketched in 
Figure U The results from the sum-ansatz (solid line) are compared with the results of the DS-improved 
superposition scheme (dashed line). 
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V. IMPROVING SELFDUALITY BY APPLYING THE ADHM FORMALISM 

The ADHM formalism [42|, has been originally formulated by Atiyah, Hitchin, Drinfeld and 
Manin in 1978 to provide a proper framework for the construction of arbitrary (anti-)selfdual solutions 
of the Yang-Mills equation of motion in flat space-time, i.e. for T = 0. The same formalism, however, 
can also be applied to obtain classical gauge field configurations at finite temperature T ^ and with 
arbitrary holonomy [J Q ■ For the convenience of the reader we will briefly review the basics of the 
ADHM formalism and then demonstrate, how the selfduality of the superposition of close calorons 
can be improved using this framework. Finally, we will again compare the improved solutions with 
the corresponding results of the naive sum-ansatz. 



A. Short summary of the ADHM formalism 



Formally, there are two construction frameworks, the Sp(N) ADHM formalism based on quaternions 
and the SU(N) ADHM formalism based on complex spinors, and since SU(2) = Sp(l) both are 
applicable for SU(2) instantons or calorons. However, it could be shown, at least for topological 
charge Q = 1, that both approaches are identical [42|,|52|]. Here, we will restrict ourselves to the Sp(l) 
construction scheme for charge Q = k > instantons or calorons. Anti-selfdual gauge fields can then 
be obtained from the corresponding selfdual gauge fields by Eq. 

The basic ingredients of the Sp(l) ADHM formalism for the construction of a charge k G N instanton 
are a k x k symmetric, invertible, quaternionic matrix B and a quaternionic row vector A. Using the 
Pauli matrices to define Oi = — iTj the space-time position x — (x, t) is encoded as a quaternion 
x = XiOi + tl allowing to define the space-time position dependent matrix B{x) by 

B{x) = B — x ■ l fe ,fc , (27) 

where \k,k is the k x k quaternionic unit matrix. The vector potential A^x) corresponding to the 
input data B and A of this construction scheme is then obtained from the quaternionic solution vector 
u{x) of the following linear system of equations 

B\x)-u{x) =A f , (28) 

which has to be solved for every point x in space-time. The vector potential is then given by the 
simple equation 

iA^x) = w\x)d^w(x) , (29) 

where a projection to the hyperspace spanned by {<j\, cr 3 } is implicit and w{x) is just the solution 
vector u{x) scaled by the potential 4>{x) according to 

w(x) = * /i • u(x) , (30) 
V <P 

<f>(x) = l+u\x) -u(x) 61. (31) 

The vector potentials obtained from Eq. (f2"9"]) solve the equation of motion in flat 1R 4 and yield the 
topological charge k if and only if the ADHM constraint 

At (x) ■ A(x) G Mat{k, k,M) , (32) 

AW =( B t>) (33> 

is fulfilled, where Mat(k, fc,IR) denotes the set of real-valued k x k matrices. In that case the k 
eigenvalues of B can be regarded as the positions of k instantons in that solution. 

One can easily see that this constraint is always fulfilled for k — 1. Therefore, all initial parameters, 
i.e. B\\ and Ai, can be chosen freely, resulting in 8 free parameters for a selfdual solution in flat 1R 4 
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with topological charge k = 1. This is the single instanton solution with B\\ defining its position and 
Ai describing its size and global SU (2) rotation. 

The time-periodic caloron, defined on IR 3 x S 1 , with topological charge Q — 1 and arbitrary 
holonomy can also be constructed by means of the ADHM formalism [l], Q. Applying the inverse 
of the gauge transformation ([3]) transforms the time-periodic caloron into a non-periodic, selfdual 
solution in flat IR 4 fulfilling the periodicity condition 

A«°(x + Pe t ) = Too ■ A^(x) ■ Vt . (34) 

It is the latter which can be constructed by the ADHM formalism. 

For the full parameter space of a single caloron the ADHM data B mn and A m , where m,n 6 Z, 
have been elaborated by Kraan and van Baal [j], Q • The corresponding ADHM data are 

/ >x r sin (2-7r(m — n)u) . . . fn _ s 

B mn = (m + ■ S m , n + p&ffi — — • (1 - <W (35) 

and 

A m = e 2mm ^C , (36) 

where £, £ are quaternions and (3 = 1 for convenience. The 8 free parameters of a charge Q = 1 caloron 
are encoded in B mn and A m . The quaternion £ describes the 4D space-time position of the caloron. 
The set of eigenvalues of B is {£ + j(3 : j G Z} and corresponds to a chain of instantons. The size 
parameter p, the spatial orientation of the monopolc axis, and the global SU(2) rotation parameter 
are determined by 

t = p-q, qeSU(2), ar = q^Tq, (37) 

where a is the direction unit vector in three dimensions between the two caloron constituents. 

This formalism can also be adopted for completely anti-selfdual instanton or caloron gauge fields, 
which can be related to their selfdual counterparts by Eq. (JB)) . The latter transformation can be built 
into the formalism by replacing Eq. (I27|) with 

B anti (x) = B - zt • 1 M . (38) 

In that case one has 

u anti {x) = u( X i), (t) anU (x) = 0(xt), and w antt (x) = u;(x f ) (39) 

leading to an exactly anti-selfdual gauge field. Its action density matches that of the corresponding 
selfdual counterpart if one inverts all spatial coordinates, i.e. x <-> . In particular, the caloron 
center £ becomes ^ for the anticaloron. 



B. Approximate multi-caloron solutions from ADHM formalism 

The ADHM formalism as described in the previous Section can in principle also be applied to 
multi-caloron configurations (all selfdual) with higher topological charge Q > 1. Thus, for a given set 
of Q desired caloron center positions £j (i — 1, Q) one is left with the problem of (i) solving the 
non-linear ADHM constraint Eq. (|32[) . (ii) fulfilling the periodicity condition Eq. (|3"4"|). and (Hi) fixing 
the eigenvalues of B according to a set of Q 'instanton chains', i.e. {eigenvalues of B} = {£; + j(3e t : 
i = 1,...,Q; j G Wi}. For some special cases, e.g. all constituent positions located along one axis, 
corresponding solutions for the ADHM data have been presented 0, 0]. But in all generality the 
ADHM constraint could not be solved due to its non-linear nature. Therefore, our idea is to construct 
approximate, multi-caloron ADHM data B PA , X PA by using the well-known single caloron ADHM 
data £?W, A^ as building blocks. Thereby, we keep insisting on the periodicity condition ii) and on 
the correct localization of the caloron positions iii) in the sense given above, while we allow for some 
violation of the ADHM constraint. However, since the ADHM data \® of each single caloron 
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solve the ADHM constraint exactly, the hope is that the newly composed multi-caloron ADHM data 
B PA , X PA can serve as an acceptable approximation for the true multi-caloron system. Later we will 
see that this pseudo-ADHM ansatz at least leads to gauge held configurations with higher selfduality 
than one would have obtained from the naive superposition of caloron gauge helds. 

Here, we obtain the pseudo-ADHM data by simply merging the ADHM data of Q single calorons 
according to 

H o • else ^ ^ m ' neZ ' M = 1,->Q, (40) 
Ka A Q+l = A$ with m G 7L, i = 1, Q . (41) 

This arrangement of the components of the single ADHM matrices Bm,n guarantees the periodicity. 
Since the set of eigenvalues of B PA is exactly the union of the eigenvalues of all flW, the caloron 
positions are unchanged by this ansatz. Additionally, the ADHM constraint still holds for every pair 
of columns of A(x) which have been filled with the ADHM data of the same caloron. On the other 
hand, for two columns belonging to different data sets i, j € {1, ...,Q} (i ^ j) one has 

(At.^(») • A PA ( X )) mQ+inQ+ . = m ,n e Z. (42) 

Therefore, the extent to which the ADHM constraint is violated is determined only by the vectors 
AW. 

This ansatz yields exact multi-caloron solutions in the case of trivial holonomy (Voo = 1) and equal 
global SU(2) rotations of all involved calorons, i.e. if all are chosen proportional to each other, 
since then one always has Xtn^Xn G H» The sizes and 4D-positions can then even be chosen arbi- 
trarily. For instantons this observation is well-known as 't Hooft ansatz (unpublished, see Ref. (53l|). 
Here, it will serve as our motivation to apply this ansatz also in the general case, where the holonomy 
is non-trivial. 

We will now see that an analytical expression for the vector potentials can be derived from this 
ansatz. Let vM'(x) denote the solution of Eq. ([28| of the i th caloron. The corresponding solution for 
B PA , X PA can then be composed as 

u%Q+i(x)=u$(x), me%, i = l,...,Q. (43) 

From Eq. (j3"T|) one derives 

Q 

<t> PA {x) = l + « t ' PA (aj)u PA (a:) = l + 53« t 'W(a;)«W(a;) 
Q 

= l + ^^W^-l) . (44) 
i=i 

The vector potentials corresponding to B PA , X PA are now obtained from Eqs. (|29[) and p0|) yielding 



Q 

E^wfw. ( 45 ) 



4> PA {x) ^ 

where A^\x) denotes the vector potential of a single caloron. Here, the other superscript indexing 
the Lie algebra components is implicit. We will refer to Eq. (|45p as the pseudo-ADHM superposition 
scheme in the following. In words, the final result of this superposition ansatz is a weighted sum of 
the single caloron gauge helds with the weighing factor </> (?) / (t> PA (x) analytically given by Eqs. (|44l 
and the analytical expression for <f>(x) which are given in [l|, |2j]. 
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Furthermore, for the real modeling of the Yang-Mills vacuum, we need the extension of this ansatz 
to compositions of calorons and anticalorons. Certainly, such configurations can neither be selfdual 
nor anti-selfdual. Nevertheless, we try such an extension. We incorporate such compositions into this 
construction scheme by replacing Eq. (127[) with 



b pa (x) = \ B ™Q+^q+j ■- xS m , n Si,j : if the object i is a caloron 

' BmQ+i,nQ+j ~ x ^m,n$i,j '■ if the object i is an anticaloron 



where m, n € 7L and i,j — l,...,Q. For a selfdual caloron the single ADHM data SW are chosen 
by Eq. (|35[) with £ = where £j encodes the desired caloron center position. For an anticaloron, 
however, we set £ = Q in Eq. (|55|) in order to locate the anticaloron at £j. 

The solution vector u PA (x) of the linear Eqs. (|28|) can then again be composed out of the single 
caloron solution vectors u PA (x). Using Eq. ([55]) we obtain 

^pa / \ _ I u m {\i,^i,x) : if the object i is a caloron 
,n Q +l I u m (Xi,^j,x^) : if the object i is an anticaloron ' 

where m G 7L and i = 1, ...,Q. Here u(X,^,x) is the solution vector of Eq. (|28|) for a single caloron 
with the center position encoded as £. Since one has 4>^(£i + x) = (f>^(£i + x^) as can be seen from 
the analytical expression of 4>(x) given in [1, 2], the potential <f> PA (x) for the extended ansatz is again 
given by Eq. (I44p . where the functions always denote the caloron ^-potential for the caloron center 
position Finally, we arrive at 



A PA (x) = 1 V<ft w (s) x ( : if the object i is a caloron ^ 

p Pj4 (a:) ^ 1 A'^ ltl ^{x) : if the object i is an anticaloron 



This is in slightly modified form the result that was already obtained in Eq. (j43)l . The only difference 

^ an 



is the appearance of the anti-selfdual vector potentials A a J ltl ^ l \x) corresponding to anticalorons at 



some center positions 

A few further remarks are in order here for the orientation of the reader. First of all, it should be 
noted that the pseudo-ADHM superposition scheme recovers the sum-ansatz (jHJ far apart from all 
caloron centers, since one has (p % '(x) — > 1 for |£j — x\ — > oo. Hence, all weighing factors become unity 
in that limit. Very close to the monopole centers, however, one has <jy^{x) 3> 1 and the corresponding 
gauge field becomes the dominant contribution to the weighted sum. More precisely, the corresponding- 
weighing factor becomes unity while all others are suppressed. Therefore, this ansatz leads to better 
selfduality within the caloron cores and recovers the sum-ansatz far away from the cores. Applying 
this ansatz to compositions of calorons and anticalorons enforces the existence of regions where the 
solution is either nearly selfdual or nearly anti-selfdual, together with regions where the field is far 
from being either selfdual or anti-selfdual. 

Finally, we demonstrate that this ansatz is a generalization of Shuryak's ratio-ansatz [5J], a tech- 
nique to improve the superposition of an instanton (x) and an antiinstanton A^ (x). Given such 
a pair with positions £12 an d instanton size parameters one directly finds 

5 (1) =6, B< a >=& \\™\=Px* (49) 



Id - A i« x \ |£i )2 - 

The separate (anti)instanton gauge field [U,[55[ is related to c^ 1 ' 2 ) by 



(50) 



A^{x) = ^Otf)l v dAog^ x \x) and A^(x) = ^Of 7^ log ^ (*) (51) 

where fj^ v , 77" denote the 't Hooft tensors and r a are the Pauli matrices. The matrices describe 
the color orientation which was encoded in A' 1:2 ). Inserting this into the pseudo-ADHM superposition 
scheme (|48|) one immediately arrives at 

a pa, , ^-tQi b v>lyivl\yi\ 4 + T f_ Q?vlA„/\y2\ 4 
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where yi t 2 = x — £12 denote the 4-dimensional, real-valued distance vectors from the (anti)instanton 
centers. This result exactly recovers the ratio-ansatz [54j j for instanton-antiinstanton pairs. 

C. Comparing pseudo-ADHM-improved solutions with the sum-ansatz 

The pseudo-ADHM superposition scheme (|48|) is now compared with the sum-ansatz (|5J) . Again, the 
continuum gauge fields A^ A (x) and AV tm (x), respectively, will be discretized on a lattice to make the 
constructed gauge fields accessible to cooling procedures. Both the improved as well as the unimproved 
configurations are then exposed to a number of cooling steps. From the cooled configurations we assess 
the strength of the alteration during the cooling process and in this way indirectly the quality of the 
applied superposition scheme. 



FIG. 8: Sketch of the caloron composition to be used for the comparison of the superposition schemes. The 
four calorons are well-dissociated (p = f.O, /3 = 1), and carry maximally non-trivial holonomy |uJ| = 0.25. 
The signs (±) correspond to the magnetic charge of the monopoles and the straight lines symbolize the Dirac 
strings. The latter are aligned along the edges of a square. We will consider the arrangement where the 
separation of two neighboring equal-charge monopoles is put equal to Az = 0.28, which is small compared to 
the 3-dimensional action lump size pzd ~ /3 = 1 of the undisturbed calorons. 

For the purpose of comparison a simple planar composition of four calorons as sketched in Figure [5] 
is considered. The calorons are well-dissociated (p = f.O, j3 = f ) and the monopole axes are aligned 
along the edges of a square such that 4 pairs of monopoles with equal magnetic charge are formed. The 
distance Az = 0.28 between two neighboring monopoles is small compared to the 3-dimensional radius 
P3D ~ /3 = 1 of their action lumps in the undisturbed case. One therefore expects strong violations 
of selfduality at all monopole positions. The alert reader might have noticed that this example is 
constructed in such a way to avoid the problem (dealt with in Section IIV[) of interaction between the 
Dirac string of one caloron with the monopole constituents of another caloron. 

Figure M shows the action density, the topological charge density, the Polyakov loop field £(x) and 
the local violation of selfduality 6(x) for the sum-ansatz. As expected, the equation of motion is 
strongly violated at the locations of the monopole pairs. This is shown by the local violation of 
selfduality S(x) or by the fact that topological charge density and action density do not match in 
these regions. From the ratio 7 of action and topological charge taken from Table HTT1 we furthermore 
see that the interaction of the close monopoles generates an action surplus of about 52%. 

In order to reduce the violation of selfduality e-cooling with e = is applied to this starting 
configuration. The corresponding cooling history is presented in Table HTT1 One sees that a relatively 
small number O(10) of cooling steps is sufficient to remove the action surplus. 

The same composition of calorons, this time generated by means of the pseudo-ADHM superposition 
scheme, is shown in Figure [TUJ Clearly, the local violation of selfduality S(x) is much smaller than 
for the previously discussed sum-ansatz. Furthermore, the topological charge and action density now 
resemble each other and the topological charge is now Q = 4 with high precision. Besides that, the 
improved selfduality can also be recognized by the smaller action surplus, being approximately 21 % 
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Action density Topological charge density 




Self-duality violation Polyakov loop field 




FIG. 9: Action density in instanton units, topological charge density, Polyakov loop field l(x) and the local 
violation of selfduality 8(x) for the gauge field configuration obtained by the sum-ansatz from the four-caloron 
composition with Az — 0.28 as sketched in Figure [8] 



Cooling steps 


Action S[A] 


Topological charge Q[A] 


1 = S/Q 





6.009 


3.967 


1.515 


1 


5.148 


3.980 


1.293 


5 


4.270 


3.985 


1.071 


10 


4.090 


3.985 


1.026 


20 


4.008 


3.984 


1.006 


50 


3.984 


3.984 


1.000 


100 


3.982 


3.984 


1.000 



TABLE III: Cooling history of the gauge field configuration obtained by the sum-ansatz from the caloron 
composition with Az = 0.28 as sketched in Figure [8] The total action S[A], the total topological charge Q[A] 
as well as the ratio 7 = S[.A]/Q[A] are shown after certain numbers of e-cooling steps with e = 0. 

as obtained from the observable 7 in Table HVl The corresponding cooling history is also presented in 
this Table. Three further remarks are in order here. 

(i) It should be noted that the action lumps obtained from the pseudo-ADHM technique are much 
sharper localized than they are for an undisturbed caloron with the same parameters. This effect 
becomes even much stronger for further decreasing monopole distances Az. The shrinkage of the 
monopoles is due to the weighing factor in Eq. (|45p with which each single gauge field is multiplied. 
As discussed in the previous Section this weighing factor guarantees that close to the monopole centers 
the corresponding gauge field dominates the superposition such that the selfduality is improved. 

(ii) The pseudo-ADHM technique leaves the monopole center positions unaltered. 

(Hi) The outcome of the pseudo-ADHM technique can not be imitated by the sum-ansatz and some 
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adequate adjustment of the parameters of the involved calorons. 



Action density Topological charge density 




FIG. 10: Action density in instanton units, topological charge density, Polyakov loop field £(x) and the local 
deviation from selfduality S(x) for the gauge field configuration obtained by the DS-improved superposition 
scheme from the caloron composition with Az — 0.28 as sketched in Figure [5] 
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4.000 
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4.260 




4.000 


1.065 


20 


4.122 




3.999 


1.031 


50 


4.022 




3.998 


1.006 


100 


4.001 




3.997 


1.001 



TABLE IV: Cooling history of the gauge field configuration obtained by the DS-improved superposition 
scheme from the caloron composition with Az — 0.28 as sketched in Figure[8] The total action S[A], the total 
topological charge Q[A\ as well as the ratio 7 = 5*^4] are shown after certain numbers of e-cooling steps 
with e = 0. 

Next, it will be demonstrated that cooling alters the gauge field configuration obtained from the 
sum-ansatz strongly, whereas the improved gauge field given by the pseudo-ADHM approach is only 
modified gently thanks to its almost perfect selfduality. Therefore, the action densities corresponding 
to both cases after applying 100 cooling steps are shown in Figure fTTl (The topological charge density 
is the same after cooling.) While the action density provided by the sum-ansatz has completely 
changed its appearance, the action density of the pseudo-ADHM improved configuration is only slightly 
altered. In the latter case the monopoles tend to drift apart from each other, which implies that the 
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monopole cores grow in their spatial extension and their maximal action density decreases. One 
may therefore conclude, that the pseudo-ADHM improved gauge field is already close to a classical 
solution by construction, under the particular condition to consist of separable monopoles. For further 
improvement (cooling) it reacts by slightly separating the monopoles which by itself are stabilized by 
the weighted superposition. 




FIG. 11: Action densities in instanton units Si for the caloron composition sketched in Figure [8] after applying 
100 e-cooling steps with e = 0. The gauge field configurations before cooling were constructed by the sum- 
ansatz (left) and by the pseudo-ADHM technique (right). 

We should draw particular attention to the special volcano-like shape of the cooled configurations 
that have developed from the sum-ansatz. Such a structure (with local minima of action density in the 
center) is characteristic for m = 2 monopoles and m = — 2 antimonopoles and has also been observed 
in analytical solutions for charge Q = 2 calorons with two equally charged monopoles coming close 
to each other while the other two were sent to infinity 6]. In contrast, the cooled pseudo-ADHM 
improved configurations do not exhibit such kind of structure. Hence, it remains unclear whether 
the pseudo-ADHM improved solutions in our case can actually be regarded as being closer to the 
exact Q = 4 solution in the configuration space, although their selfduality is certainly improved by 
far compared to the sum-ansatz. 
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FIG. 12: Dependence of the total action S[A] in instanton units Si (left) and of the total topological charge 
Q[A] (right) on the monopole separation Az for the caloron composition sketched in Figured] The results from 
the sum-ansatz (solid line) are compared to the results of the pseudo-ADHM superposition scheme (dashed 
line) . 
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Finally, the dependence of the created action surplus on the monopole separation Az is of interest. 
Figure [12] shows the total action and topological charge as measured in the pseudo- ADHM improved 
and unimproved (sum-ansatz) configurations versus Az. Since the two superposition schemes become 
equivalent at large caloron separations as discussed in the previous Section, the results for the total 
action converge to each other for Az —> oo. For Az//3 < 0.8 the action of the sum-ansatz exceeds 
the action of the improved construction. At closer separations the action surplus generated by the 
sum-ansatz becomes much larger than for the pseudo- ADHM improved configurations. The improved 
superposition scheme hence becomes meaningful at small distances, here at Az < 0.7f3. The topologi- 
cal charge is almost perfectly Q — 4, independently of the distance Az, for the improved superposition, 
and relatively constant (within a few percent) for the sum-ansatz. 

D. Comparison with an exact Q — 2 solution 

In the last section the gauge field configurations constructed by the pseudo- ADHM superposition 
scheme were shown to suffer much less from deviations from selfduality than those obtained by the 
sum-ansatz. Now we address the question whether these improved configurations are moreover also 
capable of qualitatively imitating exact higher charge solutions for some appropriately chosen single- 
caloron parameter sets. 

For that purpose we selected an exact Q = 2 solution which was provided to us by the authors 
of Q and which is shown in Fig. I13h .d.g. We then fixed the optimal single-caloron parameter sets for 
both superposition schemes independently by minimizing the square norm of the difference between the 
action densities of the exact and the superposed configurations. The optimal configuration constructed 
by the pseudo- ADHM technique is shown in Fig. [T5b.c.h. while Fig. [T5b.f.i shows the "optimal" 
configurations obtained from the sum-ansatz. 

Comparing Fig. [T3k . b, and c one clearly sees that it is indeed possible for the pseudo-ADHM 
ansatz to mimic the 4-monopole structure of the exact Q = 2 solution despite their small separations, 
while this is not the case for the sum-ansatz. For the latter ansatz the solution is completely spoilt by 
strong Dirac string interactions, which do not occur for the pseudo-ADHM technique. One should also 
notice that these Dirac string interactions can not be suppressed by the DS superposition scheme due 
to the weak caloron dissociation. Besides the strong improvements concerning the local self-duality 
also the overall measure 7 — S[A]/Q[A] for self-duality is again obviously better for the pseudo- 
ADHM ansatz. For the unimproved configuration we find 7 = 1.221 and 7 = 1.066 for the improved 
configuration. Since the self-duality aspects have already been discussed in the previous section, they 
are not explicitly demonstrated in detail here again. 

However, for the pseudo-ADHM ansatz it can clearly be observed in Fig.[T5k.h that those monopoles 
which were initially enclosed in the same caloron are still stronger connected (see Fig.[l3]h) than those 
initially belonging to different calorons (see Fig. [T5b). This asymmetry is due to the fact that there is 
no mechanism included within the pseudo-ADHM ansatz to unify the single calorons into one single, 
higher charge configuration. Therefore, this scheme will always preserve the single-caloron character 
of the superposition. This is because the single ADHM-data of the involved calorons were merged in 
the simplest possible manner. However it might be surprising to see how well this improved technique 
also works out with respect to imitating this exact Q — 2 solution in contrast to the sum-ansatz. 

VI. SUMMARY AND CONCLUSIONS 

In the present paper we have dealt with the question how approximate multi-(anti)caloron configu- 
rations can be constructed by using exact single (anti)caloron gauge fields as the underlying building 
blocks. We have discussed the first fundamental problem encountered when one attempts to carry 
over the well-known sum-ansatz, which can successfully be applied for the construction of not too 
dense multi-instanton systems, to the case of KvBLL calorons with non-trivial holonomy. The tech- 
nical answer to this problem has already been used in our recent study Q of the caloron gas where 
we have found that it works relatively well even without further improvement. This certainly applies 
only to the parameter range that we found to be relevant. In the present paper we have proposed two 
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FIG. 13: a,d,g: Action density of exact Q — 2 solution. b,e,h: Action density of pseudo-ADHM improved 
configuration. c,f,i: Action density of unimproved sum-ansatz configuration. d,e,f: Horizontal view in direction 
(dxi — l,dx3 — 1). g,h,i: Horizontal view in direction (dxi — — l,dxs = 1). 



improved superposition schemes to deal with difficulties arising when the caloron system is either too 
dense or when the calorons increasingly appear in a well-dissociated state. 

The first fundamental problem, arising from the non-vanishing A<± component due to the non-trivial 
holonomy of the considered calorons, can be overcome by superposing the caloron gauge fields in the 
algebraic gauge. This ansatz, however, implies the restriction that only calorons (and anticalorons) 
with identical holonomy become superposed, which has significant consequences when one is calcu- 
lating observables like the Polyakov loop correlator in a semiclassical, caloron based simulation. This 
is due to the long-ranged correlations induced by this superposition scheme. Whether this restriction 
has a real physical background, meaning that multi-caloron configurations with variant single-caloron 
holonomies are actually suppressed, or whether it is simply a remnant of the imperfectness of the 
applied superposition scheme, remains unanswered here. 

Secondly, we have dealt with the possible caloron-Dirac string interaction, which can easily create 
a tremendous action surplus when an additional caloron is inserted into the system in such a way 
that one of its monopoles comes too close to the Dirac string of another, well-dissociated caloron. We 
therefore have introduced the so-called DS-improved superposition scheme, which - loosely speaking - 
displaces the Dirac string before the vector potentials are superposed and moves the Dirac string back 
to the original place afterwards. By comparison with the sum-ansatz we observed a greatly improved 
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sclfduality of the configurations obtained from this (Dirac-string improved) superposition scheme. 
The excess action is essentially reduced. This already cures some, increasingly unpleasant features 
one might encounter by employing the caloron gas model in a slightly modified physical regime with 
strongly separated monopoles. 

We have checked that the DS-improved configurations are practically stable under improved cooling 
(e = cooling) in contrast to the unimproved gauge fields. From these observations we conclude that 
even without cooling the configurations constructed within the DS-improved superposition scheme 
are very close to classical configurations. This technique is also applicable in systems consisting of 
calorons and anticalorons which are not required to be close to any classical configuration and where 
cooling would be problematic as a tool to improve, in a second step, the configuration obtained by the 
sum-ansatz. The only technical restriction for the DS-improvement scheme is that the constituents 
of the dissociated caloron containing the Dirac string have to be sufficiently far separated from each 
other. We therefore conclude that the DS-improved superposition scheme might be the method of 
choice for the confinement phase if it turns out that the calorons should be more dissociated into 
monopoles (dyons) than we have found appropriate in our preceding paper Q • 

We then presented another superposition scheme suitable to improve the classicality and selfduality 
or anti-selfduality of (anti)caloron superpositions in the case of overlapping calorons or monopole cores 
(of different calorons). Our so-called pseudo-ADHM ansatz is based on the ADHM formalism. Its 
main idea is to approximate the underlying ADHM data of the actual (anti-)selfdual multi-caloron 
configuration by combining single (anti) caloron ADHM data rather than attempting to improve the 
corresponding gauge field directly. Here, we simply unify the ADHM data of all involved calorons 
or anticalorons, which guarantees (i) that the caloron center positions (subject to sampling) remain 
unchanged, (ii) that the periodicity in time is exactly preserved, and (Hi) that an analytical expression 
in terms of the single (anti)caloron vector potentials including an easily calculable weighing factor 
exists. The price to pay is to abandon to solve the (nonlinear) ADHM constraint. For trivial holonomy 
and equal orientations in color space of all calorons being involved, this ansatz yields exactly selfdual 
configurations, boiling down to the well-known 't Hooft ansatz in the case of instantons. On the other 
hand, our pseudo-ADHM superposition scheme corresponds to Shuryak's ratio-ansatz when applied 
to instanton-antiinstanton pairs. Presenting a comparison of this improved superposition scheme with 
the sum-ansatz we have demonstrated that the selfduality of gauge fields created by the pseudo-ADHM 
construction is significantly improved and the excess action essentially reduced. This becomes a major 
advantage at small caloron or monopole distances. 

The overall structure of the improved configurations as represented by the action density is relatively 
stable under improved cooling. Overlapping monopoles are driven by cooling to slightly separate from 
each other, keeping their identity as monopoles. In contrast to these observations, the configurations 
obtained from the sum-ansatz are strongly changing their structure during the cooling process. This 
means that these gauge field configurations are rather far from classical ones. An interesting property 
of unimproved (sum-ansatz) configurations is that they are not protected against the merger of equal- 
charge monopoles into higher-charge monopoles under cooling. This is witnessed by the emergence of 
a typical volcano-like structure in the action density which is known to be a characteristic feature of 
exact charge Q = 2 caloron solutions in the limit of approaching equal-charge monopoles. It should 
be kept in mind that the pseudo-ADHM ansatz is intentionally constructed in such a way to preserve 
the identity of the single monopoles. 

As a final remark, the combination of the benefits of the DS-improvement and of the pseudo- 
ADHM technique would result in a much wider applicable version of the caloron gas model that we 
have started to develop in Ref. 0]. However, since the weighing factor of the pseudo-ADHM ansatz 
changes the structure of the Dirac string, this has not been accomplished so far. This remains an 
important subject for further studies. 
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